Superconductor-proximity effect in hybrid structures: Practality versus Chaos 
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We study the proximity effect of a superconductor to a normal system with fractal spectrum. 
We find that there is no gap in the excitation spectrum, even in the case where the underlying 
classical dynamics of the normal system is chaotic. An analytical expression for the distribution of 
the smallest excitation eigenvalue E\ of the hybrid structure is obtained. On small scales it decays 
algebraically as V(E\) ~ E^ D °, where Do is the fractal dimension of the spectrum of the normal 
system. Our theoretical predictions are verified by numerical calculations performed for various 
models. 

PACS numbers: 74.45.+C, 05.45.Mt, 73.23.-b 
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During the last years there was considerable inter- 
est for the study of the statistical properties of hybrid 
normal-superconductor structures. One of the main out- 
comes of these studies was the prediction that a normal 
system in the proximity of a superconductor acquires 
characteristics that are typical of the superconducting 
state. Specifically, if the underlying classical dynamics of 
the normal system is chaotic, an energy gap in the quasi- 
particle density of states emerges above the Fermi energy. 
The mean value (Ei } of the gap was found 0, E3l to be 
proportional to the Thouless energy, 



(Ex) oc D T = GS/Att = MT5/4ir 



(1) 



where (• • •} indicates an ensemble average, M is the num- 
ber of transverse modes of the point contact between the 
superconductor and the normal system, T is the tun- 
nel probability per mode and 8 ~ 1/ L d is the mean level 
spacing of the d- dimensional normal system of linear size 
L. The product G = MT is the point contact conduc- 
tance in units of 2e 2 /h. As a matter of fact, in a recent 
publication 0] the statistical fluctuations of the lowest 
excited state around the mean- field value (E\) were stud- 
ied in the framework of the Random Matrix Theory. It 
was found that the gap distribution is a universal func- 
tion of the rescaled energy x = (E\ — (Ei))/A g , in a 
broad range |x| <C M 2 / 3 . The width of the gap distri- 
bution A g ~ M 1 / 3 ^ is parametrically smaller than the 
gap size (E\) but bigger than the mean level spacing 8 of 
the normal system. In contrast, normal systems with in- 
tegrable classical dynamics do not possess any gap near 
the Fermi energy. Instead, their density of states van- 
ishes linearly near the Fermi level. Thus, it was natu- 
rally proposed, that the appearance or not of a gap in 
the excitation spectrum of a normal system in the prox- 
imity with a superconductor can be used in the studies 
of quantum chaos as a measure for distinguishing classi- 
cally chaotic systems from integrable ones. Apart from 
the two extreme cases discussed above, a fairly good un- 
derstanding of the proximity gap was obtained also for 



the generic case of systems with mixed classical phase 
space 0]. In this case it was found that the excitation 
gap reduces below the value of fully chaotic systems QJ. 

The investigation of the proximity gap has recently 
been extended to quantum systems in the diffusive regime 
pj where it was found that a similar type of gap appears. 
The value of E\ is given by Eq. (JTJ provided that we 
substitute the ballistic conductance MT with the appro- 
priate expression DL 2 ~ d for diffusive systems. Here D is 
the classical diffusion coefficient. 

Despite all this activity, a significant class of systems 
was left unexplored. Namely hybrid structures whose 
normal part has fractal spectra. The latter exhibit en- 
ergy level statistics that are in strong contrast to the level 
repulsion predicted by Random Matrix Theory (RMT) 
ja| . Their level spacing distribution follows inverse power 
laws P(s) ~ s~P which is a signature of level clustering. 
The power (3 was found to be related with the fractal di- 
mension of the spectrum Do as (3 = 1 + Dq [7| . Realiza- 
tions of this class are quasi-periodic systems with metal- 
insulator transition at some critical value of the on-site 
potential like the Harper model El , Fibonacci chains 
0- > or quantum systems with chaotic classical limit as 
the Kicked Harper Model 

IH IHII3- Apart of their own 
interest the analysis of these systems can be illuminating 
for the understanding of the behavior of high dimensional 
disordered systems at the metal-insulator transition like 
the id Anderson model [13j ■ 

Here, for the first time, we present consequences of the 
fractal nature of the spectrum of the normal system in 
the excitation spectrum of the hybrid structure. We con- 
sider the normal system connected to the superconductor 
via point contacts supporting M channels and show that 
there is no gap in the excitation spectrum even in the case 
where the corresponding classical phase space is chaotic 
(like in the case of the Kicked Harper model). We de- 
rive an analytical expression for the distribution of the 
minimum excitation eigenvalue V(Ei), and show that its 
behavior is dictated by the fractal dimension Dq of the 
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spectrum of the normal part. Thus the nature of the clas- 
sical dynamics becomes totally irrelevant for these type of 
systems. Specifically we show that V{E{) generated over 
different Fermi energies behaves as 



1)(E~ D <> 1) 



(2) 



where E\ = 2—^ and s ma x is the maximum level spacing 
of the normal system within the energy interval that is 
used in order to generate statistics. Eq. @ is the main 
outcome of our investigation. A consequent result is that 
the mean E\ is given by 



Eh = 



1 - Da 

2- Dn 



(3) 



One has to note the lack of any dependence on the sys- 
tem size L and the number of channels M in contrast to 
Eq. Q ■ Our theoretical results are confirmed by nu- 
merical calculations performed for various systems with 
fractal spectra. 

Let us start our analysis with the Kicked Harper (KH) 
model. The system is defined by the time depended 
Hamiltonian 



H = Q 



cos(p) 



Kcos(9)^2S(t 



mT) 



(4) 



where p denotes the angular momentum, 9 the conjugate 
angle, while the kick period is T. Contrary to the stan- 
dard Harper model (corresponding to the limit K — » 0) 
this model for large enough K > 5 is classically chaotic. 

The quantum mechanics of this system is described by 
a time evolution operator for one period 

Q K 
U = cxp(— i— cosChp)) exp(— i— cos(x)) (5) 

h h 

where p = —id/d6 is the momentum operator and h~ is an 
effective Planck constant, which includes the frequency 
ratio of the unperturbed system and the external driv ing . 
For K = Q the quasi-energy spectrum is fractal [l(J, UJ$ 
and we always consider cases where Tl/2tt is strongly ir- 
rational. Using a recently proposed recipe |j| we can 
write down the corresponding quantum Andreev map T 
and find the quasi-energies of the excitation spectrum of 
the hybrid structure by direct diagonalization of T . In 
all cases considered bellow we have generated more than 
3000 values of E\ for statistical processing. 

Figure 1 shows V(E\) for various system sizes L and 
number of channels M. For small values of E\ the dis- 
tribution of the minimum excitation eigenvalue displays 
clearly an inverse power law. Moreover, it is independent 
of the number of channels and system size in agreement 
with Eq. (J2J • In the insets of Fig. 1 we also report our re- 
sults for the mean value of E\ for various system sizes and 
various numbers of channels. The inverse power law char- 
acter of the distribution V(E{) forces us to conclude that 




FIG. 1: The distribution of the lowest excited energy, gen- 
erated over an ensemble of different Fermi energies, for the 
KH model Various symbols correspond to different sys- 
tem sizes and numbers of channels. All data overlap with 
one another indicating that V(E\) is insensitive to the num- 
ber of channels and the system size. The dashed line is the 
theoretical prediction of Eq. @. In the insets we show the 
numerically evaluated (Ei) versus the system size L and the 
number of channels M. 



the probability to find a quasi-energy excitation smaller 
than (Ei) is high and thus there is no gap in the excita- 
tion spectrum (even in a probabilistic sense). 

The validity of Eqs. lL'1-ill can be verified in more cases 
in the Fibonacci chain model of a one dimensional quasi- 
crystal where various scaling exponents Do can be ob- 
tained. The normal system is described by the tight- 
binding Hamiltonian: 

H = J2 HV n (n\ + i\n)(n + l\ + \n + l)(n|) (6) 

n n 

where V n is the potential at site n. It only takes the two 
values +V and — V arranged in a Fibonacci sequence 0. 
It was shown that the spectrum is a Cantor set with zero 
Lebesgue- measure for all V > 0. The sample is in con- 
tact with M semi-infinite one-dimensional superconduc- 
tors which are attached in M randomly chosen sites. The 
quasi-energy spectrum of the hybrid structure is calcu- 
lated by employing the effective Hamiltonian approach 
0. We again find inverse power laws for the distribu- 
tions V(Ei) . Here the exponent depends on the poten- 
tial strength V, while Eq. @ still relates the resulting 
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FIG. 3: (a) The distribution V{E X ) for hybrid structures 
where the normal sample is a Fibonacci lattice with L — 987 
and M — 1 and various Vs. The dashed lines are the theo- 
retical predictions @; (b) The mean value (-Ei) versus the 
theoretical prediction (|3J for various V values. The straight 
line y = —0.039 + 0.993a; represents the best linear fit. 



FIG. 2: The distribution V(Ei) of a hybrid structure consist- 
ing of a a Fibonacci sample of size L (normal part) attached 
with M semi-infinite one-dimensional superconductors. The 
data for various L, M fall one on top of the other, indicating 
that V(Ei) is independent of L,M. In the insets we report 
the mean value (Ei) of the distribution as a function of system 
size L and number of channels M. 



statistics to the fractal dimension D E . 

In Fig. 2 we report our results for V(E±) for V = 1.4 
and various system sizes L and number of channels M . 
In the insets we also plot the mean value (Ex) of the 
distribution. Similarly with Fig. 1 we observe that V(Ei) 
and consequently (Ex) are independent of L and M. In 
Fig. 3 we summarize our results for various V values. 
A nice agreement between our numerical data and the 
theoretical predictions (|2I3|) is observed [l4| . 

The above results call for a theoretical explanation. 
Our starting point is the observation that the distribu- 
tion V(Ei) rescaled in appropriate way is the same for the 
normal and the hybrid structure. This assumption is ver- 
ified in Fig. 4 where we plot the distribution of the min- 
imum excitation level E\ for a representative case where 
the normal part is a Fibonacci lattice with V — 1.4. The 
overlap between the resulting distributions of the hybrid 
structure and of the corresponding normal system is evi- 
dent. In order to understand this phenomenon one should 
recall that the eigenstates of the normal system are not 
extended like typical chaotic eigenstates, but they have 
fractal structure. As a result, most of the eigenstates 
have intensities at the boundary with a superconductor 
which are so small that one can consider that they are 



practically not affected by the proximity of the system to 
the superconductor. Thus these eigenstates and the cor- 
responding eigenenergies are solutions of the eigenvalue 
problem for the hybrid structure as well. We point here 
that the perturbative assum ptio n used in our argument, 
was verified numerically in |l5|. where it was observed 
that the effect of the coupling of systems with fractal 
spectra to external leads is a small perturbation for the 
most of the eigenstates. 

The distribution V(E\) for the normal system can be 
derived in the following way. We scan the spectrum of 
the normal system with a set of Fermi energies which 
have resolution given by e. Then for fixed E\ we count 
the number of Fermi energies that are in a distance E\ 
with tolerance e from the next larger eigenenergy of the 
normal system. This is given by the number of level 
spacings which are larger than E\ i.e. J E ™ ax p(s)ds where 
•Smax is the maximum level spacing and p(s) is the level 
spacing distribution. The normalized gap density in the 
limit of e — > is then given by 



V(E X 



Ie7** p^ ds 



Jo Smax dEx J E ™ X p(s)ds 



(7) 



Substituting in the above equation the expression 
p(s) = ,s~( 1+I?0 ) we eventually get 

V(Ex) = (-±--l)(Ex D °-l) 



for 



(8) 



where Ex = is the rescaled energy gap. Notice that 
an additional factor 2 should be introduced in the defini- 
tion of the rescaled energy gap i.e. Ex = 2-^- once we 
turn to the distribution of the hybrid structure, due to 
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FIG. 4: The distribution V(E\) for a hybrid structure where 
the normal part is a Fibonacci system with L = 987 and 
M — 50 and V = 1.4 (o). Overplotted (*) is the corre- 
sponding distribution V(Ei) for the same Fibonacci system 
disconnected from a superconductor. 



the fact that the excitation levels for the hybrid structure 
are coming in pairs. In Figs. 1 and 3a we plot with dashed 
lines the above theoretical prediction. Using Eq. JSJ the 
mean E\ and can be easily evaluated leading to Eq. 

In conclusion we have studied the statistical proper- 
ties of the excitation spectrum of a normal system with 
fractal spectra in the proximity to a superconductor. We 
have found that the underlying classical dynamics is ir- 
relevant and that the statistical properties of the quasi- 
spectrum depends only on the fractal nature of the nor- 
mal system. Such a system can be realized by a two- 
dimensional electron gas subject to a perpendicular mag- 
netic field and periodic potential in contact with a su- 
perconductor |16( . It is known that the normal system of 
this type can be mapped onto the Harper model |7llslll4j 
possessing fractal spectrum at the critical point. 

Finally, our results can be of interest for the quan- 
tum optics community with respect to studies of re- 
flection of light by a dielectric medium in front of a 
phase-conjugation mirror |l7^ . This problem is the op- 
tics analogue of Andreev reflection A wave inci- 
dent at frequency loq + Alu is retro-reflected at frequency 
uj — Alu where cu is the pumped frequency of the phase- 
conjugation mirror. The analogue of luq and of the fre- 
quency shifts Aw are the Fermi energy Ep and the excita- 
tion energies E respectively. The latter were shown here 
to be controlled by the fractal dimension Do. Keeping in 
mind the above analogies it is reasonable to conjecture 
that Alu is controlled as well by Dq in quasi-periodic op- 
tical structures like e.g. in Fibonacci quasicrystals [19(. 
Thus controlling Dq we can tune Alu which was predicted 
in 0| to affect drastically the reflected intensity from a 
phase conjugation mirror for a certain parameter range. 
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